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Over  the  last  60  years,  a  considerable  amount  of  time  and  money  has  been  spent 
improving  our  understanding  of  combustion  instability  in  solid  rocket  propulsion  systems. 
Over  this  period,  significant  knowledge  has  been  accumulated  that  can  influence  the  acoustic 
stability  of  solid  propulsive  systems.  Unfortunately,  many  of  these  lessons  learned  about 
combustion  instability  remain  the  knowledge  of  a  select  few  in  government  and  industry  who 
work  with  combustion  instability  on  a  daily  basis.  This  paper  attempts  to  organize  many  of 
these  “rules  of  thumb”  that  propellant  formulators  and  motor  designers  need  to  be  aware  of 
in  order  to  minimize  the  chances  of  combustion  instability.  In  addition,  several  mathematical 
relationships  are  presented  which  can  be  used  to  predict  the  frequency  of  potential  acoustic 
modes  and  determine  resultant  thrust  oscillations  produced  by  acoustic  oscillations.  Also 
included  are  some  key  fundamental  equations  which  can  be  used  to  gain  insight  into 
combustion  instability  in  rocket  motor  systems  and  there  is  a  short  section  on  motor 
instrumentation.  This  paper  is  not  an  attempt  to  be  an  exhaustive  study  into  combustion 
instability,  but  rather,  will  attempt  to  list  in  a  clear  fashion  some  of  the  more  important 
lessons  learned  and  empirical  observations  of  solid  propellant  combustion  instability.  This 
paper  emphasizes  composite  propellants,  but  many  observations  apply  to  double  base  and 
composite  modified  double  base  propellants,  as  well. 
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Thermal  Diffusivity 

V 
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Acoustic  Velocity 

D 

Non-Dimensional  Frequency 

V 

= 

Volume 

I.  Introduction 

Solid  propellant  combustion  instability  is  the  amplification  or  attenuation  of  acoustic  oscillations  by  solid 
propellant  combustion  processes  in  a  rocket  motor.  Combustion  instability  in  solid  rocket  propellant  motors  has 
been  a  continuing  problem  since  the  first  rockets  were  used  in  World  War  II.  Shortly  after  the  war,  research  was 
begun  trying  to  understand  the  phenomenon. 1Pncc  Pressure  oscillations  can  take  the  form  of  longitudinal,  tangential 
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Longitudinal  Oscillations 


Tangential  Oscillations 


Radial  Oscillations 


and,  more  rarely,  radial  oscillations.  See  Figure  1.  The 
term  linear  instability  refers  to  oscillations  which  are 
sinusoidal  in  nature  and  can  be  linearly  decomposed 
into  discrete  sin  waves.  The  term  non-linear  instability 
refers  to  oscillations  which  contain  many  acoustic 
modes  and  are  often  characterized  by  steep  fronted  non- 
sinusoidal  waveforms.  These  oscillations  cannot  be 
linearized  into  discrete  sin  waves.  This  type  of 
instability  is  often  caused  by  injecta  or  debris  passing 
through  the  nozzle  which  pulse  the  motor.  Non-linear 
oscillations  are  often  accompanied  by  DC  pressure 
shifts  and  limiting  amplitudes  of  the  oscillations. 


Figure  1.  Solid  Rocket  Motor  Acoustic  Mode  Types. 


Normally,  combustion  instability  pressure 
oscillations  are  on  the  order  of  a  few  percent  of  the 
mean  pressure  and  cause  no  problems  with  rocket  motor  performance.  Occasionally,  however,  these  oscillations  can 
reach  higher  levels  which  can  cause  DC  pressure  shifts  in  the  mean  pressure,  cause  coupling  effects  with  guidance 
and  control  sections,  structural  coupling,  or,  in  the  worst  case,  catastrophic  motor  failure.  An  example  of  a  motor 
experiencing  severe  oscillations  accom¬ 
panied  by  large  DC  pressure  shifts  and 
subsequent  motor  failure  is  shown  in  Figure 
2.  Dozens  of  solid  rocket  motors  developed 
over  the  last  60  years  have  experienced 
some  form  of  combustion  instability.2'3 
These  motors  experienced  everything  from 
minor  linear  acoustic  oscillations  to  violent 
non-linear  oscillations  accompanied  by 
large  ballistic  pressure  shifts  and  in  some 
cases  motor  destruction.  In  many  of  these 
systems  considerable  cost  was  added  to  the 
development  or  the  programs  were 
canceled.  In  others,  after  careful 
considerations  of  the  consequences,  it  was 
decided  to  live  with  non-destructive 
oscillations. 


Combustion  Instability. 


During  development,  solid  rocket 
motors  can  experience  different  types  and 
degrees  of  combustion  instability.  Due  to 
these  problems,  considerable  research  has  been  performed  during  the  last  60  years  to  investigate  combustion 
instability.4  Theories  were  developed  in  attempts  to  model  primary  physical  mechanisms  which  influence 
combustion  instability.  These  theories  were  at  first  purely  one-dimensional  in  nature,  but  more  recently  have 
branched  off  into  multi-dimensional  and  non-linear  models  of  combustion  instability.  A  brief  overview  of  one¬ 
dimensional  models  will  be  presented  in  order  to  reinforce  some  of  the  empirical  trends  observed  and  provide  the 
reader  with  some  background.  Empirical  rules  of  thumb  to  prevent  or  at  least  minimize  effects  of  combustion 
instability  have  been  developed  through  experiences  with  many  different  solid  motor  systems.  Some  of  these  will  be 
presented  in  this  paper.  Other  techniques  have  also  been  developed  over  the  years  to  predict  at  what  frequencies 
instability  might  occur  and  what  are  the  thrust  oscillations  induced  by  pressure  oscillations.  Many  of  these 
relationships  will  also  be  presented. 


Hopefully,  the  content  of  this  report  will  be  of  use  to  the  propulsion  community.  An  effort  was  made  to  reference 
many  of  the  observations  reported  herein,  however,  many  important  references  may  have  been  inadvertently  omitted 
which  can  reinforce  (or  add  doubt)  to  the  conclusions.  I  encourage  all  interested  persons  to  report  to  me 
discrepancies,  missing  observations  and  missing  references.  Please  send  your  comments  or  suggestions  to  Dr.  Fred 
S.  Blomshield,  Code  474200D,  M/S  6204,  Naval  Air  Warfare  Center,  China  Lake,  CA  93555-6106,  760-939-3650, 
fred.blomshield  @  navy.mil. 
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II.  Linear  Stability  Prediction 


A.  Driving  Mechanisms 


1)  Pressure-Coupled  Response'.  The  pressure-coupled  response  is  the  amplification  or  coupling  of  the  combustion 
processes  taking  place  at  the  propellant  surface  with  acoustic  pressure.  It  is  the  often  referred  parameter  to  describe 
combustion  instability  characteristics  of  a  propellant.  The  pressure-coupled  response  is  mathematically  defined  as 
the  ratio  between  the  perturbed  burning  rate  over  the  mean  burning  rate  to  the  perturbed  pressure  over  the  mean 
chamber  pressure: 


R 


pc 


rl  r 

p/p 


(1) 


The  non-dimensional  Rpc  is  a  function  of  burning  rate,  pressure,  frequency  and  propellant  physical  and 
mechanical  properties.  Zero  frequency  corresponds  to  the  steady  state  condition  and  the  pressure-coupled  response 
is  the  pressure  exponent  assuming  the  propellant  burning  rate  obeys  the  classical  St.  Roberts  burning  rate  law.3 

r  =  cpn  (2) 

Typical  response  curves  taken  by  the  T-burner  are  shown  in  Figure  3. 6  These  curves  all  have  typical  shapes  for 
composite  propellants.  This  particular  plot  compares  four  similar  propellants  only  differing  by  three  percent  of  one 
ingredient.  At  zero  frequency  Rpc  is  the  exponent  and  as  frequency  goes  up  so  does  the  response  to  a  maximum 
value  of  around  1.25  to  2.5  depending  upon  the  propellant.  As  frequency  gets  even  higher,  the  response  falls  off  to 

zero.  Double  base  propellants  show  a  similar 
behavior  except  the  fall-off  to  zero  usually  occurs 
at  a  much  higher  frequency.  The  T-burner  is  the 
most  common  way  to  measure  response  but 
numerous  other  devices  can  be  used.7,8  Other 
devices  to  measure  the  pressure-coupled  response 
include  the  rotating  valve,  magnetic  flow  meter, 
microwave  technique  and  the  radiant  heat  flux 
device. 9’10,11’12’13 


There  are  several  reasons  for  obtaining  the 
pressure-coupled  response  of  a  propellant.  First,  it 
can  be  used  to  compare  propellants  destined  for 
the  same  application.  Once  expected  frequencies 
are  known  for  a  solid  motor  (this  will  be 
discussed  later),  the  propellant  with  the  lowest 
response  will  be  less  likely  to  drive  instabilities, 
other  factors  remaining  constant.  It  is  very 
important  to  remember  that  combustion  stability  in  a  rocket  motor  is  a  system  dependent  phenomenon  that  depends 
on  many  factors  such  as  pressure,  geometry,  structure,  nozzle  type,  etc.  The  second  reason  for  obtaining  the 
response  is  that  it  can  be  used  by  motor  stability  prediction  programs  to  compute  the  net  driving  by  propellant 
combustion  in  a  rocket  motor. 


Figure  3  Typical  Pressure- Coupled  Response  Curves 


2)  Velocity-Coupled  Response:  The  velocity-coupled  response  is  the  amplification  or  attenuation  of  combustion 
processes  taking  place  at  the  propellant  surface  by  the  acoustic  velocity.14  Unlike  the  pressure-coupled  response 
which  has  broad  support  for  its  basic  definition  and  nature,  the  velocity-coupled  response  definition  does  not  share 
this  universal  support.  Most  experts  do  agree  that  velocity  coupling  does  correspond  to  the  non-steady  component  of 
erosive  burning.  Recent  research  indicates  that  velocity  coupling  is  a  non-linear  phenomenon  and  that  conventional 
linear  approaches  cannot  be  applied.4'13'16 17  Unlike  acoustic  pressure,  acoustic  velocity  not  only  varies  in  the 
longitudinal  direction  but  also  in  the  radial  direction  due  to  a  thick  acoustic  boundary  layer.  Because  of  these 
variations,  the  acoustic  velocity  can  have  phase  differences  in  both  radial  and  longitudinal  directions.  How  the 
mechanisms  of  these  ideas  this  can  be  applied  to  motor  stability  is  beyond  the  scope  of  this  report.  Fortunately, 
velocity-coupled  response  driving  is  believed  to  be  very  small  for  many  motors.  However,  some  combustion 
instability  problems,  particularly  non-linear  problems,  indicate  that  velocity  coupling  can  be  an  important  driver  of 
combustion  instability. 
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3.  Distributed  Combustion :  The  third  driving  mechanism  is  called  distributed  combustion  response.  Distributed 
combustion  response  is  the  interaction  of  the  acoustic  field  with  burning  metal  particles  taking  place  away  from  the 
propellant  surface.  This  phenomenon  is  only  a  factor  when  dealing  with  highly  metallized  propellants,  and  only 
occurs  in  those  whose  burning  characteristics  allow  the  metal  to  bum  in  a  distributed  manner  throughout  a  motor 
chamber.  There  have  been  several  research  efforts  examining  this  phenomenon.18'19  Because  it  occurs  in  a  relatively 
small  percentage  of  motors,  it  will  not  be  discussed  in  any  great  detail  in  this  paper. 

B.  Damping  Mechanisms 

1.  Nozzle  Damping :  The  principle  damping  mechanisms  in  a  solid  motor  are  nozzle  damping,  particle  damping, 
mean  flow/acoustic  interactions  and  structural  damping.  Nozzle  damping  is  usually  the  largest  damping  mechanism 
in  a  motor,  particularly  with  longitudinal  and  mixed  transverse/longitudinal  modes.  Conceptually,  it  is  very  simple. 
As  acoustic  pressure  waves  come  in  contact  with  the  nozzle  throat  of  a  motor,  some  of  this  energy  is  transmitted 
through  the  throat  and  radiated  out  to  the  environment.  The  theoretical  prediction  of  nozzle  damping  is  very  well 
established.20'21’22’23 

2.  Particle  Damping:  Particle  damping  only  applies  to  propellants  containing  solid  species  in  the  combustion 
products,  such  as  metallized  propellants  or  those  with  acoustic  stability  additives.  The  amount  of  particle  damping  is 
very  dependent  on  the  mass  fraction  of  particles  in  the  flow  field  and,  most  importantly,  particle  size.  Larger 
particles  damp  lower  frequencies,  while  smaller  particles  damp  out  higher  frequencies.  The  calculation  of  particle 
damping  is  also  well  established. 5’20’24’25  Figure  4  shows  relative  particle  damping  versus  particle  size  for  six 
different  frequencies.  This  plot  was  computed  using  mono-sized  particles  of  the  same  mass  fraction.5'26'27  It  is 
important  to  note  how  higher  frequencies  are  much  more  efficiently  damped  than  lower  frequencies  for  the  same 
mass  fraction  of  particles.  Figure  5  is  a  useful  curve  that  indicates  optimum  particle  size  to  damp  a  particular 
frequency.  Both  Figure  4  and  Figure  5  are  for  typical  solid  propellants  and  are  reasonably  insensitive  to  formulation. 


Figure  4  Relative  Particle  Damping  versus  Particle  Figure  5  Optimum  Damping  Particle  Size  versus 
Size  at  Six  Frequencies.  Frequency. 

3.  Flow  Interactions :  Mean/acoustic  flow  losses  are  those  losses  associated  with  mean  flow  interactions  with  the 
acoustic  field.  There  are  generally  two  schools  of  thought  as  to  how  this  occurs  for  flow  near  the  surface  of  the  pro¬ 
pellant  in  the  boundary  layer.  The  first  is  called  flow  turning  losses.28'29  Conceptually,  this  is  the  work  that  is  done  to 
incoming  flow  as  it  acquires  acoustical  motion.20’21’23  Combustion  products  enter  normal  to  the  surface  and  must  be 
accelerated  to  the  axial  direction.  That  change  of  direction  is  the  origin  of  the  term  "flow  turning."4  The  second  is 
called  boundary  layer  losses  ’  ’  “’  It  is  also  a  form  of  analysis  to  examine  the  acoustic  boundary  with  flow  normal 
to  the  surface.  It  is  sometimes  known  as  the  "Flandro  boundary  layer  term."  Understanding  the  phenomenon  often 
deals  with  vorticity  and  acoustic  boundary  layers.  Another  type  of  flow  interaction  which  is  separate  from  the  above 
terms  is  vortex  generation  of  flow  around  sharp  corners.34'35’35'37  Large  segmented  boosters  can  have  this  type  of 
acoustic  coupling.38’39  The  final  type  of  flow  interaction  is  referred  to  as  wall  losses.  These  are  frictional  losses 
resulting  from  flow  over  exposed  case  walls.  Often  this  term  is  small,  since  usually  there  is  little  exposed  case  walls. 

4.  Structural  Damping:  The  final  damping  mechanism  is  structural  damping.  As  the  name  implies,  this  is  damping 
due  to  deformation  of  the  motor  structure  due  to  acoustic  pressure  oscillations.  The  motor  structure  can  include  the 
actual  case  as  well  as  the  motor  liner.  Normally,  this  is  thought  to  be  very  small  compared  to  other  losses  in  all  but 
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very  large  rocket  motors  where  acoustic  loading 
can  be  great.  Very  elastic  exposed  case  liners  can 
also  make  this  term  significant.  A  pictorial 
representation  of  all  of  the  gains  and  losses  in  a 
solid  rocket  is  shown  in  Figure  6. 


C.  Motor  Stability 


Linear  stability  theory  has  been  successfully 
applied  to  the  prediction  of  solid  propellant 
rocket  motors.  The  linear  assump¬ 

tion  implies  that  the  change  of  acoustic  pressure 
amplitude  in  a  motor  can  be  expressed  as  an 
exponential.  The  following  equation  can  be  used: 

P  =  pQea“t  (3) 

Here,  P  is  the  acoustic  wave  amplitude,  P0  is 
the  initial  amplitude  value,  aM  is  rate  of  acoustic 
energy  change  in  the  motor  and  t  is  time.  If  aM  is  negative,  the  motor  is  stable,  if  aM  is  positive,  the  motor  is 
unstable.  The  stability  of  a  solid  rocket  motor  is  determined  by  the  net  sum  of  acoustic  gains  and  losses  in  the 
system.  The  linear  assumption  further  assumes  that  the  contributions  to  motor  stability  can  be  simply  added  and 
subtracted  to  determine  a  motor’s  total  stability,  aM.  The  linear  relationship  for  total  motor  alpha  is: 

aM  =  aPC  +  aVC  +  aDC  +  aNOZ  +  aPART  +  aMF  +  aSD  (4) 

In  Eq.  (4),  aPc,  aVc  and  Ctix:  is  the  driving  terms  due  to  pressure  coupling,  velocity  coupling  and  distributed 
combustion,  respectively;  and  aN,)y,  aPART,  aMF  and  <%/,  are  the  damping  terms  due  to  nozzle  damping,  particle 
damping,  mean  flow  interactions  and  structural  damping,  respectively. 

Sample  stability  versus  time  plots  for  a  reduced  smoke  solid  rocket  motor  are  shown  in  Figure  7  and  Figure  8. 
Figure  7  shows  the  total  motor  stability  alpha,  at  three  different  motor  pressures  for  the  same  motor  geometry 
and  propellant.  Figure  8  shows  the  driving,  damping  and  total  stability  alphas  for  a  single  motor.  These  curves  were 
calculated  using  the  Standard  Performance  Prediction/Standard  Stability  Program  (SPP/SSP)  code  developed  under 
both  Air  Force  and  Navy  sponsorship  and  recently  improved  to  predict  non-linear  motor  behavior.47'48'49'50'51  The 
inputs  to  this  program  include  complete  motor  geometry  and  detailed  propellant  properties,  including  the  propellant 
response  function.  The  program  uses  a  grain  design  and  ballistic  code  to  provide  inputs  to  the  stability  module.  The 
program  allows  prediction  of  linear  motor  stability  as  a  function  of  time  with  a  single  set  of  inputs. 


Figure  7.  Total  Motor  Stability  at  Three  Pressures.  Figure  8.  Detailed  Motor  Stability  Map. 


Figure  6.  Acoustic  Energy  Gains  and  Losses  in  a  Solid 
Rocket  Motor. 
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D.  Stability  Observations 

The  following  section  describes  some  theoretical  observations  or  trends  for  solid  propellants  and  solid  rocket 
motor  combustion  instability  which  have  a  theoretical  background.  Many  of  the  following  statements  have 
experimental  evidence  to  reinforce  the  theory. 

Looking  at  total  motor  stability  in  Figure  7  and  Figure  8,  the  difference  between  the  value  of  motor  stability  and 
zero  damping  is  often  referred  to  as  the  limit  or  margin  of  linear  stability.  For  many  motors  the  magnitude  of  this 
margin  becomes  smaller  as  the  motor  burns.  This  implies  that  the  tendency  of  a  motor  to  become  unstable  often 
increases  with  time.36,40,41,45’46'47'52  Referring  to  Figure  8,  it  can  be  seen  that  the  magnitude  of  both  driving  alphas 
(those  that  are  positive)  and  damping  alphas  (those  that  are  negative)  decrease  in  magnitude  as  the  motor  burns  to 
completion.  This  can  be  explained  by  examining  the  following  expression: 

^(...)dS  (Usually  Constant) 

Burning  Surfaces 
aX  =  - r - 

|(...W  (5) 

Internal  Volume  (Always  Increases  During  Burn) 

Except  for  particle  damping,  all  driving  and  damping  alphas  can  be  written  in  this  form.  The  subscript  "X"  refers 
to  various  terms.  The  burning  surface  area  of  most  motors  is  reasonably  constant,  while  the  internal  volume  is 
always  increasing  as  propellant  is  consumed.  This  implies  that  as  the  motor  burns,  the  magnitude  of  the  various 
alphas,  |  a  %  |  ,  will  become  smaller.  For  metallized  systems  where  particle  damping  is  significant,  total  motor 

stability  approaches  the  particle  damping  alpha  as  the  motor  burns  to  completion.  This  is  because  the  level  of 
particle  damping  for  typical  motors  is  constant  due  to  a  constant  burning  surface  area.  As  a  general  rule,  if  a  reduced 
smoke  propellant  rocket  is  to  be  unstable,  it  will  be  toward  the  end  of  burn.  If  a  metallized  propellant  rocket  motor  is 
to  be  unstable,  it  will  be  earlier  in  burn  compared  to  a  non-metallized  system.  Figure  9  depicts  this  graphically. 


As  a  general  rule  of 
thumb,  motor  designer  should 
not  design  a  motor  with  the 
majority  of  burning  area  in  aft 

end.40, 41, 45,46,52,53  A  ^  aft 

geometry  is  an  example.  The 
driving  mechanisms  of  velo¬ 
city  coupling  and  distributed 
combustion  have  terms  like 
the  following: 

j)  pvu(...)ds  (6) 

Burning  Surface 


a 


motor 


Time 


Figure  9.  Total  Motor  Stability  Alpha  versus  Time  for  Systems  with  and 
without  Particle  Damping. 


At  the  aft  end  of  the  motor,  both  mean  velocity,  u  ,  and  oscillating  pressure,  p  ,  become  large.  ( v  is  the  oscil¬ 
lating  velocity.)  Because  of  these  terms,  the  magnitude  of  velocity  coupling  and  distributed  combustion  terms  can 
be  large  in  magnitude.  This  will  tend  to  make  a  motor  less  stable.  Nearly  all  stability  integrals  have  a  p  in  them.  For 

longitudinal  oscillations,  p  is  very  close  to  zero  near  the  center  of  the  motor  from  front  to  back.  Ideally,  most  of  the 
surface  area  should  be  at  the  center  of  the  motor  with  little  burning  surface  at  the  ends  to  minimize  driving  stability 
integrals.  Unfortunately,  this  makes  it  very  difficult  to  remove  the  casting  mandrel  using  conventional  motor  casting 
techniques.  For  these  reasons,  the  forward  end  is  the  recommended  location  for  most  of  the  burning  surface  area. 


Raising  the  mean  operating  pressure  in  a  motor  will  tend  to  make  the  motor  less  stable.  Raising  the  mean 
pressure  decreases  both  nozzle  damping  and  flow  turning  losses.  This  can  lead  to  lower  margins  of  stability  and 
small  changes  in  the  propellant  response  can  lead  to  an  unstable  motor.  Figure  7  shows  the  effect  on  total  motor 
stability  as  chamber  pressure  is  increased. 
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III.  Empirical  Observations 


The  following  observations  are  primarily  empirical  in  nature,  but  often  have  some  theoretical  basis  as  well. 
Although  they  are  grouped  by  category,  although  it  is  important  to  remember  that  some  items  could  be  listed  under 
multiple  categories. 


A.  Burning  Rate,  Pressure  and  AP  Particle  Size 

Various  observations  about  burning  rate,  pressure  and  ammonium  perchlorate  (AP)  particle  size  can  be  made  as 
to  their  effects  on  combustion  instability.  These  effects  are  grouped  together,  as  they  are  often  interrelated. 

1.  Very  fine  and  ultra  fine  (less  than  1  micron)  AP  oxidizer  particles  usually  give  high  pressure-coupled  response 
within  a  given  propellant  burning  rate  group.54  55  Propellants  with  fine  AP  have  been  known  to  cause  significant 
combustion  instability  behavior  in  motors.  Particles  less  than  5  microns  in  diameter  should  be  used  with  caution. 

2.  In  general,  if  AP  particle  size  is  constant,  a  higher  burning  rate  will  give  a  lower  pressure-coupled  response.54 
This  appears  to  be  at  odds  with  (1)  above.  If  AP  particle  size  is  held  constant,  higher  burning  rates  can  be  obtained 
by  the  use  of  different  binders,  with  additives  and/or  with  catalysts. 

3.  Very  coarse  AP,  i.e.,  greater  than  200  microns,  can  and  lead  to  motor  instability  problems.  In  addition,  there  is 
some  evidence  the  coarse  AP  causes  an  increase  in  the  non-linear  acoustic  erosivity. 56,57 

4.  The  propellant  response  function  shows  a  strong  dependence  on  mean  pressure,  although  the  pressure  at  which 
the  propellant  response  function  is  a  maximum  changes  among  different  propellant  types  and  burning  rates.54  The 
response  function  dependence  on  pressure  is  well  known,  but  at  what  pressure  the  response  is  a  maximum  at  a 
particular  frequency  can  only  be  determined  by  some  form  of  response  testing.  The  most  common  form  of  testing  is 
the  T-burner. 


5.  Propellants  with  a  low  pressure  exponent,  n  in  Eq.  (2),  will,  in  general,  have  a  lower  response  function. 
Propellants  with  zero  or  negative  pressure  exponents  such  as  mesa  burning  propellants,  are  very  desirable  from  a 
combustion  instability  point  of  view.58  The  lower  the  pressure  exponent,  the  better. 


6.  In  general,  as  mean  motor  pressure  is  increased,  other  factors  remaining  the  same,  a  motor  will  tend  to  become 


less  stable. 


40,41,45,46,52,55,59 


Figure  10. 
Instability. 


Raising  mean  motor  pressure  reduces  the  margin  of  stability  due  to  decreases  in  both 

nozzle  and  flow  turning  losses,  see  Figure 
7.  Furthermore,  a  stable  motor  will  tend  to 
be  more  susceptible  to  non-linear  pulsed 
instability.40,45,46, 52,59’60,61  An  example  of 
this  is  shown  in  Figure  10.  This  plot 
presents  the  ballistic  pressure  traces  for 
three  identical  motors  except  for  nozzle 
throat  sizes  which  modify  the  motor 
chamber  pressure.42,45  46  These  motors 
were  all  pulsed.  The  lowest  pressure 
motor  could  not  be  pulsed  into  instability. 
The  middle  pressure  motor  was  pulsed 
late  in  burn  and  non-linear  longitudinal 
instability  was  triggered  with  DC  pressure 
shifts.  The  highest  pressure  motor 
experienced  violent  spontaneous 
tangential  mode  instability  that  caused 
extreme  chamber  pressure  increases  and 
eventually  motor  failure.  All  three  motors 
had  the  same  geometry  and  propellant. 


Example  of  Motor  Pressure  Effect  on  Combustion 
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To  summarize  these  statements,  if  one  increases  burning  rate  with  catalysts,  propellant  combustion  response  will 
tend  to  go  down.  If  one  increases  burning  rate  with  fine  AP,  propellant  combustion  response  will  tend  to  go  up.  Very 
fine  or  very  coarse  AP  is  not  good  from  a  combustion  instability  point  of  view.  These  effects  can  be  explained  by 
the  following;  for  very  fine  AP  crystals  burning  in  a  fuel  binder,  chemical  reaction  processes  are  kinetically 
controlled  and  have  a  relatively  high  reaction  order.  The  reaction  order  is  the  pressure  exponent  in  the  chemical 
reaction  equilibrium  equation.  For  very  large  AP  crystals,  combustion  is  controlled  by  an  AP  mono-propellant  flame 
which  also  has  a  high  reaction  order.  In  between  the  extremes  of  particle  size,  combustion  processes  are  believed  to 
be  controlled  by  more  diffusional  effects  which  are  not  as  sensitive  to  pressure  oscillations.62 

B.  Flow  Effects 

Flow  effects  due  to  internal  geometry  of  a  solid  rocket  motor  can  contribute  to  the  stability  of  a  motor.  The 
following  observations  help  illustrate  this  point. 

1.  The  design  of  a  motor  should  minimize  flow  around  sharp  edges.  These  can  lead  to  vortex  shedding,  which  can 
couple  with  internal  motor  acoustics.34'35'36  Avoid  protruding  grain  inhibition  or  actual  propellant  geometries  in 
which  the  flow  is  constricted  in  the  longitudinal  direction.  Avoid  radial  slots  like  those  used  for  stress  relief. 

2.  Oscillations  created  by  vortex  shedding  are  usually  small  in  amplitude,  i.e.,  less  than  1  percent.  Vortex  shedding 
is  a  more  common  problem  in  large  segmented  boosters  like  the  space  shuttle  solid  rocket  motors.  In  these  motors, 
internal  motor  flow  around  sharp  corners  between  segments  leads  to  vortex  shedding.  Oscillations  in  the  SRM  and 
RSRM  are  believed  to  be  caused  by  coupling  between  these  large  scale  vortices  and  acoustic  modes  of  the  motor 
chamber.63'64  Because  of  their  size,  these  frequencies  are  typically  very  low,  less  than  50  Hz.  As  mentioned  above, 
these  oscillations  are  low  in  magnitude,  however,  it  must  be  remembered  that  a  single  psi  of  pressure  oscillation  can 
lead  to  hundreds  of  thousands  of  pounds  of  thrust  oscillations  in  large  motors.  How  to  determine  these  thrust 
oscillations  will  be  discussed  later 

3.  The  grain  design  should  not  have  an  abundance  of  burning  area  in  the  aft  end.40'41'45,46,52’53  This  can  increase 
driving  due  to  velocity  coupling  and  distributed  combustion  for  metallized  propellants.  This  can  also  decrease 
nozzle  damping  if  more  of  the  aft  case  is  exposed  as  the  motor  burns. 

4.  Multiple  nozzles  should  be  more  stable  due  to  increased  nozzle  damping.  However,  the  design  must  be  done 
correctly.  A  large  flat  aft  end  surface  among  the  nozzles  can  reflect  acoustic  waves  very  efficiently.56  The  area 
between  the  nozzles  must  be  designed  as  to  not  reflect  acoustic  waves.  The  nozzle  entrance  angle  is  also  an 
important  parameter.  There  is  experimental  evidence  that  a  submerged  nozzle  with  reverse  cavity  flow  slightly 
reduces  the  effective  nozzle  damping.65  Because  nozzle  damping  is  one  of  the  principle  damping  mechanisms,  small 
changes  in  aft  geometry  and  nozzle  type  and  placement  can  lead  to  big  changes  in  motor  stability. 

5.  For  star,  wagon  wheel,  or  other  similarly  slotted  geometries,  the  number  of  star  points  should  be  an  odd  number. 
An  odd  number  of  slots  will  make  it  less  likely  for  the  geometry  and  internal  flow  to  encourage  tangential  mode 
oscillations.  A  four  point  star,  for  instance,  might  encourage  a  second  tangential  mode.  However,  recent  evidence 
suggests  that  this  may  not  necessarily  be  true  and  that  the  number  of  slots  makes  no  difference.66 

6.  L*  (L-Star)  is  a  type  of  non-acoustic  instability  which  often  occurs  near  the  beginning  of  motor  burn 
immediately  after  ignition  when  the  chamber  volume  is  small  and  the  chamber  pressure  is  low.  As  its  name  implies, 
this  form  of  non-acoustic  combustion  instability  does  not  correspond  to  an  acoustic  mode  of  a  rocket  motor.  Non¬ 
acoustic  combustion  instability  results  from  a  coupling  between  the  combustion  and  the  blow-down  or  characteristic 
time  of  the  rocket  chamber.  Thus,  in  an  acoustic  wave,  the  pressure  varies  both  in  time  and  space,  through  the 
chamber,  according  to  an  acoustic  mode.  In  contrast,  during  a  manifestation  of  non-acoustic  combustion  instability, 
the  pressure  is  constant  throughout  the  chamber,  varying  only  in  time.  The  sensitivity  to  L*  instability  (sometimes 
called  chuffing  or  sputtering)  is  related  to  motor  volume  over  the  nozzle  cross  sectional  area,  L*  =  V/AT.67  The 
lower  the  L*,  the  more  likely  a  motor  will  experience  L*  instability.  Figure  11  is  a  compilation  of  many  L*  values 
for  a  variety  of  propellants  and  indicates  both  stable  and  unstable  combustion  regions.57  This  figure  is  a  general 
guide  to  when  a  propellant  in  a  motor  might  be  unstable. 
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C.  Particle  Effects 


Particles  in  the  combustion  products,  both 
inert  and  burning,  can  have  a  significant 
effect  on  combustion  instability.  The  effect  of 
particles  on  acoustics  has  been  studied  for 
many  years,  and  its  specific  application  to 
solid  rockets  has  evolved  over  the  past  forty 
years.  Inert  particles  include  aluminum  oxide, 
AI2O3,  produced  from  aluminum  combustion 
and  stability  additives  such  as  zirconium 
carbide,  ZrC,  or  inert  A1203  added  to  the 
propellant  during  mixing.  These  particles 
cause  the  particle  damping  mentioned  above. 
Stability  additives  are  crucial  ingredients  in 
nearly  every  modern  reduced  smoke, 
composite  propellant  formulation  in  use 
today.  They  have  been  proven  very  effective 
in  reducing  and/or  eliminating  combustion 
instability  in  many  solid  rocket  systems.  In  one  recent  study  numerous  reduced  smoke  composite  propellant  tactical 
sized  rocket  motors  were  intentionally  pulsed. 6,405:  Many  of  the  motors  containing  propellants  without  additives 
were  pulsed  into  violent  unstable  non-linear  instability,  accompanied  by  extensive  DC  pressure  shifts,  and  often 
resulted  in  nozzle  or  motor  case  failure.  These  same  motors,  cast  with  propellants  containing  one  percent  ZrC,  8|im 
AI2O3  or  90pm  AI2O3  in  place  of  one  percent  the  ammonium  perchlorate  (AP),  could  not  be  pulsed  into  unstable 
behavior  even  when  pulsed  twice  as  hard.  Figure  12  shows  this  amazing  effect  graphically.  Burning  metal  particles 
can  cause  a  distributed  combustion  coupling  with  internal  motor  acoustics  and  drive  instabilities.  For  most 
metallized  propellants,  metal  particles  burn  very  close  to  the  propellant  surface  and  do  not  have  time  to  couple  with 
the  acoustics.  However,  evidence  suggests  that  these  burning  metal  particles  can  contribute  to  combustion  instability 
for  high  burning  rate  propellants  at  high  pressure  through  the  mechanism  of  distributed  combustion.  Below  are  some 
observations  about  combustion  instability  behavior  of  particles  in  solid  rocket  motors. 


Figure  11.  Non-Acoustic  Combustion  Instability  Boundary 
for  a  Variety  of  Solid  Propellants. 


1.  Particle  damping  at  a  particular 
frequency  is  very  dependent  on  particle 
size  (see  Figure  5). 


2000 


6,26,27 


2.  Given  equal  mass  fractions  of 
particles,  smaller  particles  are  much  more 
efficient  at  damping  higher  frequencies 
than  larger  particles  are  at  damping  lower 
frequencies  (see  Figure  4). 6 ,26'27 

3.  Some  additives  are  effective  in 
reducing  combustion  instability  behavior 
of  solid  propellants.  This  is  done  by  the 
reduction  of  the  propellant  response 
function  and  by  increasing  particle 
damping. 


6,40,40,52,53,60 


Large  Amplitude 
Non-Linear 
Oscillations 


Motor  Failure 

Motor  without 
Stability 
Additive 

Identical  Motor 
I  %  ZrC 


4.  Particulates  in  the  flow  are  very  non¬ 
linear  in  nature. 40’52’60,68  This  means  that 
linear  theory  does  not  predict  the 


2  3  4 

Time  -  seconds 

Figure  12.  Motor  Stability  with  and  without  Stability  Additive. 

observed  behavior  of  particles  in  solid  rocket  motors.  The  very  nature  of  particle  damping  causes  one  particle  to 
damp  out  many  different  acoustic  modes,  adding  to  the  non  linearities. 


5.  The  effect  of  distributed  combustion  on  high  burning  rate,  metallized  propellants  at  high  pressure  may  not  solely 
be  due  to  an  increase  in  acoustic  driving  by  burning  aluminum  particles,  but  also  a  lowering  of  particulate  damping 
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at  the  desired  frequency  due  to  a  reduction  in  the  normally  damping  A1203  particles.  These  burning  particles  are 
damping  as  well,  but  at  different  frequencies.  The  net  distributed  combustion  result  is  less  particle  damping  at  the 
desired  frequency  of  oscillation.18 

6.  For  high  burning  rate  propellants  at  high  pressure,  there  is  a  direct  relationship  between  aluminum  particle  size 
and  combustion  response.  The  smaller  the  aluminum  particle  size,  the  higher  the  response.18  For  conventional 
burning  rate  propellants  at  lower  pressure,  there  is  some  evidence  that  smaller  aluminum  particles  may  yield  a  lower 
response.55’56  There  is  also  evidence  that  the  type  of  aluminum,  i.e.,  flake  or  particulate,  as  well  as  aluminum 
coatings,  such  as  grease,  can  effect  the  combustion  response  of  propellants.55'56  To  be  certain,  aluminum  particle 
size  and  type  effects  on  combustion  response  must  be  measured  for  different  classes  of  propellants. 

7.  Metallized  propellants  produce  a  bimodal  or  trimodal  distribution  of  A1203  particles  that  will  be  present  in  the 
combustion  chamber.18  The  smaller  size  fraction  is  often  around  1  micron  and  is  directly  produced  by  the  burning 
aluminum  particles  in  the  propellant.  The  larger  size  fraction  is  produced  by  metal  agglomerations  at  the  surface 
and  is  often  as  large  as  200  microns.  The  acoustic  damping  is  strongly  dependent  upon  the  particle  size,  see  Figure 
4  and  Figure  5.  The  relative  mass  fraction  of  the  fine  and  coarse  size  fractions  is  usually  determined  experimentally 
by  burning  a  propellant  sample  in  a  particle  collector.  After  careful  size  analysis  the  particle  damping  can  be 
computed. 

D.  Chemistry  Effects 

This  area  is  very  important  in  determining  combustion  instability  behavior  of  propellants.  Unfortunately,  it  is 
also  the  area  that  lacks  good  experimental  evidence  to  understand  the  controlling  mechanisms  that  cause  some 
propellants  to  respond  to  pressure  oscillations  and  some  not  to  respond  to  pressure  oscillations.  The  following 
observations  are  based  mainly  on  theoretical  predictive  modeling  and  are  somewhat  speculative. 

1.  A  propellant  with  an  energetic  pressure  dependent  diffusion  flame  is  more  likely  to  have  instability  problems. 
AP  propellants  typically  have  this  form  of  flame.  CL- 20  and  ADN  oxidizer  propellants  have  this  type  of  flame  to  a 
lesser  extent.  HMX  oxidizer  propellants  do  not  have  this  type  of  flame  and  would  typically  drive  less  instabilities 
than  AP  based  propellants.69 

2.  Preliminary  results  indicate  that  CL-20  and  ADN  will  have  very  weak  diffusion  flame  with  inert  binders.  This 
will  result  in  very  small  particle  size  effects  on  burning  rate  (like  HMX).70  This  probably  indicates  that  CL-20  and 
HMX  particle  size  variations  will  have  little  impact  in  combustion  instability  behavior.69 

3.  Initial  evidence  shows  that  ADN  has  a  low  melting  point,  low  surface  temperature  and  low  pressure  exponent 
which  tends  to  indicate  a  significant  condense  phase  reaction  on  the  surface.  In  general,  the  greater  the  condensed 
phase  reactions,  the  higher  the  temperature  sensitivity  which,  in  turn,  could  lead  to  instability  behavior.57 

4.  The  higher  the  temperature  sensitivity,  the  more  likely  a  propellant  will  be  subject  to  combustion  instabi¬ 
lity.  1,Ku  ota  Of  the  observations  in  this  section,  this  one  has  the  most  substantiated  experimental  and  theoretical  basis. 


E.  Non-Linear  Effects 

The  area  of  non-linear  instability  is  an  area  in  which  recent  progress  has  been  made.  There  have  been  numerous 
papers  on  the  subject  covering  both  experimental  observations  and  modeling.4’30,31’32’33'40,41’45’46’49’50’51’52’53’59’60’61’66’68’ 

’  ’  The  details  of  this  work  is  well  beyond  the  scope  of  this  report.  Some  important  conclusions  and  comments 
about  non-linear  combustion  instability  can  be  made.  Non-linear  behavior  is  a  function  of  hardware  (motor  or  test) 
with  dependence  on  burning  rate  and  response  function  values  over  all  frequencies.  The  term  non-linear  implies 
many  different  harmonic  modes  are  all  excited  simultaneously  and  the  system  cannot  be  linearized  into  discrete 
wave  forms.  This  type  of  instability  is  often  caused  by  injecta  or  debris  passing  through  the  nozzle  which  pulse  the 
motor.  Non-linear  oscillations  are  often  accompanied  by  DC  pressure  shifts  and  oscillation  limiting  amplitudes. 

1.  Non-linear  stability  is  dependent  upon  propellant  response  over  the  entire  frequency  range.4’15’40'45'46’52’68  Non¬ 
linear  combustion  instability  is  believed  to  be  related  to  energy  feedback  from  higher  harmonics  to  lower  harmonics. 
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In  Figure  13  propellant  “A”  has  a  higher 
response  peak  while  propellant  “B”  has  a  much 
broader  response  peak  that  stays  high  as 
frequency  increases.  Propellant  “A”  should  be 
more  stable  from  a  non-linear  point  of  view 
since  propellant  “B”  will  be  more  susceptible 
to  amplifying  non-linear  high  frequency  waves. 

2.  Pulsing  a  motor  can  give  quantitative 
information  to  determine  how  stable  that  motor 
is.  Pulsing  can  be  used  to  determine  the  margin 
of  stability.  If  this  knowledge  is  known,  motor 
developers  can  feel  much  more  confident  about 
changing,  for  example,  AP  particle  size  or 
manufacturer  without  causing  the  motor  to 
become  unstable.40'45'46'52  Good  examples  of 
motor  pulsing  are  shown  in  Figure  2  and 
Figure  12. 

3.  Stability  additives  can  be  very  effective  in 
eliminating  pulsed  non-linear  instability.  6’40,52’60’68  Figure  12  is  an  example  of  this.  The  additives  work  by  reducing 
the  combustion  response  over  a  broad  frequency  range  and  by  increasing  the  particle  damping.  Figure  13,  although 
only  an  illustration,  shows  how  a  stability  additive  might  reduce  the  response  at  higher  frequencies.  In  this  example 
propellant  “A”  would  have  the  stability  additive.  The  pulse  amplitude  required  to  trigger  non-linear  instability  will 
be  increased  and  the  DC  shift  and  instability  wave  amplitude  decreased  with  the  addition  of  stability 
additives.40'45’46'52'60'61 


Figure  13.  Illustration  of  Response  Function  effect  on  Non- 
Linear  Stability. 


4.  One  parameter  that  has  been  used  to  rank  propellants  in  laboratory  devices  as  to  their  susceptibility  to  non-linear 
instability  is  as  follows'74'75'76 


n  = 


(7) 


In  this  expression,  a  is  the  linear  growth  rate  of  an  oscillation,  f  is  the  frequency,  R x  is  the  limiting  amplitude 
of  an  oscillation  and  p  is  the  mean  pressure. 


5.  The  susceptibility  to  pulsed  instability  increases  with  increasing  chamber  pressure.40’45'46'52'60'61  Figure  10 
illustrates  this.  Furthermore,  the  triggering  pulse  amplitude,  i.e.,  the  pulse  amplitude  needed  to  cause  instable 
combustion,  decreases  as  the  motor  pressure  is  increased.45'52  59  68 


6.  In  general,  the  lower  the  motor  mean  gas  velocity,  the  more  likely  a  motor  can  be  pulsed  or  triggered  into  non¬ 
linear  instability.15’17’40’45’46’52'77  Lower  mean  velocities  result  from  higher  motor  operating  pressures  and  happen  later 
in  burn  due  to  a  greater  cross-sectional  area  of  the  chamber. 

7.  The  amount  of  flow  reversal  in  the  forward  end  of  a  motor  from  a  pulse  can  contribute  to  the  triggering  behavior 
of  a  pulse  to  cause  non-linear  instability.57  61  This  is  related  to  (6)  above,  since  more  flow  reversal  can  occur  when 
the  mean  flow  velocity  is  lower. 


8.  There  is  direct  relationship  between  DC  pressure  shifts  and  limiting  amplitude  of  the  oscillations  for  a  motor 
experiencing  pulsed  non-linear  instability.40'45'46'52  Figure  14  supports  this  assumption.  It  is  important  to  note  that  the 
slope  of  the  line  in  Figure  14  is  dependent  upon  propellant,  geometry,  chamber  pressure  and  other  motor  specifics. 

9.  Oscillatory  decay  rates  from  motor  pulses  are  independent  of  pulse  amplitude40'45’46’52’59’68’77  This  is  only  true  up 
to  the  triggering  pulse  amplitude  at  which  time  the  motor  will  go  unstable.  The  growth  rate  of  unstable  oscillations 
will  then  remain  independent  of  pulse  amplitude  for  pulses  greater  than  the  triggering  amplitude.  High  pulse 
amplitudes  can  excite  instability  in  an  otherwise  stable  motor. 
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IV.  Commonly  Performed  Computations 

Often  when  determining  the  character  of  an 
instability  occurring  in  a  motor,  simple  calculations 
are  performed  to  help  determine  what  kind  of 
instability  it  is.  Other  common  calculations  are 
performed  to  determine  effects  of  combustion 
instability  on  motor  performance.  Finally,  some 
calculations  are  performed  which  seem  rather 
unusual,  but  are  performed  often  enough  to  be 
reported  here.  Also,  a  section  on  the  T-Burner  and 
how  it  works  is  included. 


A.  Frequency  Analysis 


Figure  14.  Relationship  between  DC  Pressure  Shifts 
and  Limiting  Amplitude. 


With  a  very  rudimentary  knowledge  of  a  motor’s 
interior  geometry,  various  natural  acoustic 


frequencies  can  be  computed.  From  this  information  and  a  frequency  analysis  of  an  oscillation,  it  is  usually  possible 
to  determine  how  the  motor  is  oscillating. 


1.  Longitudinal  Modes:  The  most  common  form  of  instability  in  motors  with  composite  propellants  results  in 
longitudinal  acoustic  waves  traveling  from  the  head  end  of  a  solid  motor  to  the  aft  end  and  back  again,  see  Figure  1 . 
These  oscillations  always  have  pressure  anti-nodes  at  both  ends  of  the  motor.  For  the  first  acoustic  mode,  the 
pressure  node,  i.e.,  where  no  oscillations  are  occurring,  is  at  the  center.  The  second  mode  has  an  additional  anti-node 
at  the  center  and  nodes  at  the  1/4  and  3/4  point  along  the  longitudinal  axis  of  the  motor.  Location  of  the  pressure 
nodes  and  anti-nodes,  except  the  end  anti-nodes,  will  vary  somewhat  depending  upon  internal  geometry  of  the 
motor.  The  location  of  acoustic  velocity  nodes  and  anti-nodes  is  the  exact  opposite  of  pressure  nodes  and  anti-nodes. 
To  approximate  the  frequency  of  longitudinal  modes,  the  following  equation  may  be  used: 


Where:  n 
a 

L 


Acoustic  mode  number,  i.e.,  1  for  first  mode,  2  for  second  mode,  etc. 

Actual  speed  of  sound  in  motor  cavity,  approximate  values  are  2900  and  3500  ft/sec  for 
metallized  and  reduced/minimum  smoke  propellants,  respectively 
Chamber  length  in  feet 


2.  Transverse  modes:  Transverse  acoustic  modes  include  both  tangential  and  radial  mode  oscillations,  see  Figure  1. 
Tangential  mode  oscillations  are  most  common  in  double  base  propellants  and  hardly  ever  seen  in  metallized 
propellants.  Radial  modes  are  very  seldom  seen  in  solid  rockets  and  more  commonly  seen  in  liquid  propulsion 
rockets  engines.  Combination  modes  are  also  quite  rare,  but  do  happen  occasionally.  Nodes  and  anti-nodes  for 
tangential  modes  are  not  well  defined,  since  these  modes  can  actually  rotate  with  in  the  motor.  Internal  geometry 
such  as  slots,  however,  can  cause  these  modes  to  stay  stationary.  A  frequency  spectrum  of  a  transverse  mode  will 
show  a  decrease  in  frequency  as  the  motor  bums,  since  internal  diameters  are  getting  larger.  The  following  equation 


can  be  used  to  approximate  transverse  modes: 

( ®  p ,  t  ) a 


Table  1.  Tangential  and  Radial  Mode 
Coefficients 


2  R 


(9) 


Where: 


a 

R 

a 

T 

P 


PS 


Speed  of  sound,  see  Eq.  (8) 

Internal  radius  at  time  of  instability 
mode  coefficient,  see  Table  1 

Tangential  mode  number 
Radial  mode  number 


p 

X 

0 

1 

2 

3 

0 

0 

1.220 

2.233 

3.238 

1 

0.586 

1.697 

2.714 

3.726 

2 

0.972 

2.135 

3.173 

4.192 

3 

1.337 

2.551 

3.611 

4.643 

Example:  Determine  first  tangential  mode  coefficient,  p-  0  and  T=  1,  then  the  coefficient  to  use  is,  etoj=0.586. 
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3.  Omega,  Non  Dimensional  Frequency :  The  non-dimensional  frequency,  omega-Fi,  is  often  used  by  theoreticians 
to  understand  the  nature  of  the  instability.1'78'79  Many  composite  propellants  whose  pressure-coupled  response  has 
been  evaluated  by  the  T-burner  have  a  response  peak  with  an  omega  value  of  between  5  and  30  and  often  has  a 
value  around  10. 57  80  It  is  the  ratio  of  the  acoustic  time  to  the  thermal  conduction  time.  It  is  computed  by  the 
following: 


^  _  l7rfv 

~  -2 
r 

Where:  f  frequency  in  cycles  per  second 
r)  thermal  diffusivity 

r  burning  rate  in  units  of  in/sec 


(10) 


In  this  expression,  T)/r2  is  the  characteristic  time  of  the 
thermal  wave.  Some  typical  values  of  thermal  diffusivity  in 
units  of  in2/sec  can  be  found  in  Table  2. 


Table  2.  Typical  Thermal  Diffusivity  Values. 


Propellant  Family 

Value 

PBAN  AP 

2.64  x  10-4 

Polyurethane-no  A1 

1.82  x  10-4 

Pure  AP 

3.45  x  10-4 

PBAN  based  with  10%  A1 

3.25  x  10-4 

Aluminized  Polyurethane 

2.80  x  10-4 

Aluminized  HTPB  or  PBAN 

1.60  x  10-4 

B.  Thrust  Oscillations  Due  To  Pressure  Oscillations 


Once  a  motor  does  experience  pressure  oscillations,  it  is  often  desirable  to  know  what  impact  these  oscillations 
will  have  on  thrust  of  the  motor.  Many  people  incorrectly  assume  that  all  one  needs  to  do  is  multiply  the  nozzle  area 
by  the  peak-to-peak  pressure  oscillation.  This  is  incorrect,  and  will  yield  much  smaller  oscillations  than  actually  will 
occur.  It  is  also  incorrect  that  the  severity  of  the  thrust  oscillations  will  be  the  same  as  the  pressure  oscillations,  i.e., 
the  magnitude  of  the  oscillation  over  the  mean  value.  In  fact,  thrust  oscillation  severity  is  almost  always  worse  than 
that  of  the  pressure  oscillations.  The  following  two  approximate  methods  assume  longitudinal  oscillations. 
Transverse  modes  usually  cause  no  thrust  oscillations.  The  first  method  is  a  quick  and  dirty  method  which  is  always 
very  conservative.  It  will  yield  larger  thrust  oscillations  than  actually  occur.  It  provides  a  good  upper  limit  to  thrust 
oscillations.  It  is  as  follows: 

A F  =  2  Ac  AP  (11) 

Where:  AF  Thrust  oscillations 

Ac  Chamber  cross-sectional  area 
AP  Peak  to  peak  pressure  oscillations 


The  factor  of  "2"  arises  from  the  fact  that  pressure  oscillations  act  not  only  on  the  nozzle  but  on  the  aft  and 
forward  end  of  the  motor  as  well.  Since  pressure  anti-nodes  at  the  motor  ends  are  180°  apart,  a  positive  delta 
pressure  acts  on  one  end  and  a  negative  delta  pressure  acts  on  the  other  end.  A  more  exact  method  is  as  follows: 

A F  =  [F  +  PAAE  -  (PN  +  PH  )AC  (12) 

‘N 


Where: 


AF  Thrust  oscillations 
AP  Peak-to-peak  pressure  oscillations 
F  Mean  thrust  at  ambient  pressure 
Ae  Nozzle  exit  area 

Ac  Chamber  cross-sectional  area  at  time  of  oscillation 
Pa  Ambient  pressure 
p^j  Nozzle  end  pressure 
Pjj  Head  end  pressure 


Equation  (12)  includes  the  following  assumptions:  frictionless  flow,  constant  cross  sectional  area  in  motor, 
oscillations  at  head  and  aft  end  equal  and  180  phase  between  head  and  aft  ends.  One  should  notice  that  as  a  motor 
bums,  thmst  oscillations  become  worse  for  the  same  pressure  oscillations  because  the  chamber  cross  sectional  area 
is  increasing.  An  alternate  form  of  this  equation  is  sometimes  written  as  thrust  change  per  psi  of  oscillation: 
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(13) 


AF_  =  F  +  PaAe  -  (PN  +  PH)AC 

A P  PN 

An  example  of  how  to  use  these  thrust  oscillation  equations  is  as  follows: 


Assume: 


F  =  10000  lbs 

Ae  =  20  in2 

A  F 

10000  lbs  +  ( 1 0  psi )  (20  in  )  -  (990  psi  +  1010  psi )  80  in  ~ 

Ac  =  80  in2 

A P 

990  psi 

PA  =  10  psi 

PN  =  990  psi 

=  151-^ 

=  1010  psi 

psi 

This  means  that  for  every  single  psi  of  pressure  oscillation  the  motor  will  produce  150  pounds  of  thrust 
oscillation.  The  short  method  would  say  160  pounds  per  oscillatory  psi. 


C.  Combustion  Gas  Velocity  Calculation 


The  combustion  gas  velocity  is  the  velocity  at  which  the  gaseous  solid  propellant  combustion  products  move 
perpendicularly  away  from  the  surface.  The  computation  of  this  parameter  is  important  to  estimate  the  velocity  of 
burning  metal  particles  leaving  the  surface  to  determine  if  distributed  combustion  is  taking  place.  To  approximate 
gas  velocity  from  a  burning  propellant,  first  an  assumption  is  made  that  the  gases  follow  the  ideal  gas  law  and  the 
density  and  the  velocity  of  the  gas  are  given  by: 


P  = 


PM 

RT 


vgas 


r propellant 


(  \ 

P solid 

yP solid  y 


(14,15) 


In  general,  the  approximate  gas  velocity  can  be  given  by  the  following  table.  The  further  you  are  from  the  base 
pressure  and  burning  rate  given  in  Table  3,  the  more  the  error  will  be  in  these  approximations.  The  injection  Mach 
number  is  typically  around  0.0029. 


Table  3.  Gas  Ejection  Velocity  Approximations 


Gas  Velocity 

Propellant  type  and  burning  rate  condition 

v  =  23858  (r/p) 

Metallized  propellant  at  1000  psi  and  0.5  in/sec,  v  =  1 1.93  ft/sec 

v  =  20075  (r/p) 

Reduced  smoke  propellant  at  1000  psi  and  0.5  in/sec,  v  =  10.04  ft/sec 

v  =  24550  (r/p) 

High  burning  rate  metallized  propellant  at  3000  psi  and  5.0  in/sec,  v  =  40.92  ft/sec 

D.  T-Burner  Theory  and  Pressure-Coupled  Response  Determination 

The  standard  way  to  measure  the  combustion  response  is  by  the  T-burner.7  The  T-Burner  has  been  around  in  one 
form  or  another  since  1958  and  is  shown  in  Figure  15.  Two  disks  of  propellant  of  equal  thickness  are  placed  on  each 
end  of  a  1.5-inch  diameter  pipe  combustor  and  they  are  ignited  simultaneously.  Ideally,  they  will  also  burn  out 
simultaneously.  The  geometry  of  the  T-Burner  was  designed  to  provide  an  environment  which  is  ideally  suited  to 
study  the  effect  of  acoustic  pressure  oscillations  on  solid  rocket  propellant.  Assuming  a  1st  longitudinal  mode,  like 
the  one  shown  in  Figure  1  is  present  in  the  burner,  the  maximum  driving  in  the  T-Burner  is  at  the  ends  where  the 
pressure  oscillations  are  a  maximum  and  all  the  propellant  is  located.  The  acoustic  velocity,  cross  flow  velocity  and 
mean  flow  are  all  zero  at  the  ends  as  to  isolate  the  pressure-coupled  response.  The  test  frequency  depends  on  the 
burner  length  and  the  combustion  gas  temperature.  Data  are  obtained  by  pulsing  the  burner  during  the  propellant 
burn  and  after  burn  out.  The  pressure  amplitude  rate  of  change  of  the  oscillations  is  measured  by  a  piezoelectric 
quartz  pressure  transducer.  To  determine  when  burnout  occurs,  phototransistors  are  located  at  each  end  to  measure 
the  light  output  at  burnout.  For  some  propellants,  pulsing  is  not  required  and  the  T-burner  acoustic  oscillations  grow 
spontaneously.  In  either  case,  the  difference  between  the  alpha  during  sample  burn,  a>,  and  the  decay  alpha  after 
burnout,  a2,  is  known  as  the  combustion  alpha,  ac. 
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High  Frequency  Pulser 

Pressure  Transducer  (at  each  end) 


Propellant  Cross  Section 
(1.5  in  Dia,  0.25  in  thick) 


T-Burner  Lengths  Vary 
From  5  to  50  Inches 
Which  Yields  Frequencies 
of  4000 -300  Hz 


Figure  15.  T-Burner  Diagram 


a.Q  =«]-«£,  while  «£>  =  a2(/l)  (16,17) 

Where:  c*i  pressure  decay  rate  constant  during  burn 
cx2(fi)  pressure  decay  constant  after  burning  with 
correction  to  frequency  of  ai 


From  the  computed,  etc,  burner  length  and  propellant 
properties,  the  combustion  response  is  computed  from: 


P  pc  ~ 


acP 


4/ apprb{SB/Sc) 


(18) 


Where:  P 
f 

P_p 

rb 

SB/Sc 

a 

<Gi 


mean  pressure 

frequency 

propellant  density 

measured  burning  rate 

propellant  burning  surface  area  to  channel 

area  ratio 

theoretical  gas  speed  of  sound 
measured  speed  of  sound,  anl  =  2fL, 

L  =  burner  length 


V.  Motor  Instrumentation 


This  section  is  just  to  remind  motor  designers  and  test  personnel  that  combustion  instability  in  rocket  motors  can 
only  be  examined  if  the  test  motor  is  instrumented  correctly.  It  is  strongly  recommended  that  all  development 
motors  be  instrumented  to  observe  combustion  instability.  Even  if  no  instability  is  measured,  test  records  provide  a 
starting  point  if,  at  some  point  in  the  motor’s  lifetime,  a  problem  does  occur.  The  following  lists  important  points  for 
good  motor  instrumentation  for  examining  combustion  instability  behavior.  Reference  40  contains  an  excellent 
review  of  motor  instrumentation  as  well  as  data  reduction  methods. 


1.  Pressure  data  is  always 
preferred  over  strain  gage  or 
accelerometer  data. 

2.  Piezoelectric  high  freq¬ 
uency  quartz  crystal  pressure 
gages  should  be  used 
whenever  possible.  The  only 
exception  to  this  is  for  very 
large  motors,  like  the  Space 
Shuttle  SRMs  where  very  high 
frequency  response  is  not 
needed.  Figure  16  depicts  a 
suggested  motor  instrumenta¬ 
tion  data  path  for  high 
frequency  pressure  measure¬ 
ment. 


■  ) 


(^Remote  Ground^) — 


—  LP  Filter 


High  Gain 
AC  Data 


Channels  1  &  2: 

Low  Gain  AC  Data 
For  High  Level  Signals 


Channels  3  &  4: 

High  Gain  AC  Data 
For  Low  Level  Signals 


Channels  5  &  6: 
DC  Data  For 
High  Level  Signals 


Figure  16.  Motor  Pressure  Instrumentation  Data  Path 


3.  A  conventional  strain-gage 
low  frequency  pressure 
transducer  should  also  be  used  to  monitor  steady-state  chamber  pressure  and  accurately  measure  any  possible  DC 
pressure  shifts. 


4.  The  best  place  to  measure  longitudinal  modes  is  at  either  end  of  the  motor  where  the  acoustic  pressure  anti-nodes 
are  located. 
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5.  The  best  place  to  measure  transverse  modes  is  one  gage  mounted  at  0°  and  one  at  120°  around  the  motor  case 
near  where  the  oscillation  is  occurring.  This  is  often  impossible;  however,  so  gages  mounted  at  the  end  of  the  motor, 
near  the  outside  radius  at  0°  and  120°  may  be  used.  Often,  the  end-mounted  gages  will  work,  but  measured 
amplitudes  may  be  wrong. 

6.  Whenever  possible,  use  redundant  gages.  Full-scale  motor  firings  are  too  expensive  not  to  take  the  extra 
insurance. 

7.  Pressure  gages  should  be  as  close  to  the  motor  chamber  as  possible,  i.e.,  very  close-coupled  to  the  chamber 
pressure.  High  frequency  gages  mounted  on  extension  tubes  can  loose  their  frequency  response  and  the  use  of 
extension  tubes  can  add  acoustic  modes  that  are  not  part  of  or  even  related  to  actual  chamber  oscillations.  The 
resonant  frequency  of  a  cavity  installation  is  inversely  proportional  to  the  passage  length  and  can  be  determined 
from  the  following  equation.81 


4  (L  +  0.85 Z)) 

In  this  equation,  a  is  the  speed  of  sound  in  the  cavity,  L  is  the  passage  length  and  D  is  the  passage  diameter.  For 
short  cavities  that  are  on  the  order  of  the  transducer  diameter,  this  resonant  frequency  is  very  high  usually  does  not 
interfere  with  measurements  of  the  acoustic  pressure  oscillations  taking  place  in  the  motor  chamber.  As  the  cavity 
length  increases,  sharp  fronted  pressure  waves  may  excite  the  passage  to  resonance,  an  organ  pipe  effect,  and  mask 
the  true  behavior  pressure  oscillations. 


VI.  Conclusions 

This  report  has  brought  together  many  facets  of  combustion  instability:  theoretical  issues  dealing  with  driving 
and  damping  mechanisms;  motor  stability  prediction  and  various  theoretical  observations,  and  empirical 
observations  about  oxidizer  particle  size,  propellant  burning  rate,  mean  flow,  particle,  chemistry  and  non-linear 
effects  were  all  addressed.  Also  presented  were  various  commonly  performed  computations,  including  how  to 
estimate  a  motor's  longitudinal  and  transverse  acoustic  modes,  thrust  oscillation  prediction,  gas  velocity  calculations 
and  a  brief  description  on  how  a  T-burner  works.  In  addition,  a  short  section  dealing  with  motor  instrumentation  to 
monitor  combustion  instability  was  discussed.  It  is  hoped  that  this  paper  will  aid  motor  designers  and  motor  firing 
data  analysis  to  better  understand  what  causes  and,  most  importantly,  how  to  prevent  combustion  instability. 
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